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Abstract — We consider multi-antenna cooperative spectrum 
sensing in cognitive radio networks, when there may be multiple 
primary users. A noise-uncertainty-free detector that is optimal 
in the low signal to noise ratio regime is analyzed in such a 
scenario. Specifically, we derive the exact moments of the test 
statistics involved, which lead to simple and accurate analytical 
formulae for the false alarm probability and the decision thresh- 
old. Simulations are provided to examine the accuracy of the 
derived results, and to compare with other detectors in realistic 
sensing scenarios. 

Index Terms — Cognitive radio; multi-antenna spectrum sens- 
ing; multiple primary users; locally best invariant test. 

I. Introduction 

In Cognitive Radio (CR) networks, dynamic spectrum 
access is implemented to mitigate spectrum scarcity. Namely, a 
secondary (unlicensed) user is allowed to utilize the spectrum 
resources when it does not cause intolerable interference to 
the primary (licensed) users. Spectrum sensing is the first key 
step towards this dynamic spectrum access scenario. 

Prior work on cooperative spectrum sensing predominately 
employ the assumption of a single active primary user. Based 
on this assumption, several eigenvalue based sensing algo- 
rithms have been proposed recently [1-8]. These algorithms 
are non-parametric, i.e. they do not require information of the 
primary user, in contrast to e.g. feature detection. Moreover, 
they achieve optimality under different assumptions on the 
knowledge of the parameters. The assumption of a single 
primary user is made as the investigations in the literature have 
mainly focussed on CR networks, where the primary users are 
TV or DVB systems. In these systems the single active primary 
user assumption is, to some extent, justifiable. In addition, 
assuming a single primary user leads to analytically tractable 
problems. 

The single primary user assumption may fail to reflect the 
situation in forthcoming CR networks, where the primary 
system could be a cellular network, and the existence of more 
than one primary user would be the prevailing condition. Using 
existing single primary user detection algorithms in such a 
scenario will induce performance loss. Despite the need to 
understand multiple primary user detection, the results in this 
direction are rather limited. A heuristic detection algorithm 
based on the ratio of the extreme eigenvalues is investigated 



in [9, 10], but its detection performance turns out to be sub- 
optimal [6]. An optimal detection algorithm in the presence of 
multiple primary users, based on the spherical test, has been 
proposed in [11] and subsequently studied in [12]. However, 
numerical evidence suggests that this detector does not per- 
form particularly well when Signal-to-Noise Ratio (SNR) is 
relatively low [12]. Spectrum sensing in the low SNR regime 
is a practical and challenging issue in cooperative spectrum 
sensing. For example, recent FCC regulations require that the 
secondary devices must be able to detect signals with SNR as 
low as —18 dB [13]. To address this challenge, in this paper 
we consider a multiple primary user detector that is optimal in 
the low SNR regime. In particular, we investigate its detection 
performance by deriving a closed-form moment expression of 
the test statistics. Using the derived moments, approximations 
to the false alarm probability and the decision threshold are 
constructed. The derived approximations are easily computable 
and simulations show that they are accurate for the considered 
sensor sizes and number of samples. 

The rest of this paper is organized as follows. In Section II 
we propose the optimal low SNR detector for multiple primary 
user sensing after outlining the signal model. Performance 
analysis of the proposed detection algorithm is addressed in 
Section III. Section IV presents numerical examples to verify 
the derived results and to study the detection performance in 
diverse scenarios. Finally in Section V we conclude the main 
results of this paper. 

II. Problem Formulation 

A. Signal Model 

Consider the standard model for if-sensor cooperative de- 
tection in the presence of P primary users, 

x = Hs + crn (1) 

where x e C K is the received data vector. The K sensors 
may be e.g. K receive antennas in one secondary terminal 
or K secondary devices each with a single antenna, or any 
combination of these. 1 The K x P matrix H — [hi, ... , hp] 
represents the channels between the P primary users and the 

'This collaborative sensing scenario is more relevant when the K sensors 
are in one device, since for multiple collaborating devices, accurate time 
synchronization between devices are needed and communications to the fusion 
center becomes an issue. 



K sensors. The P x 1 vector s = \s\, . . . , Sp]' denotes zero 
mean transmitted signals from the primary users. The K x 
1 vector an is the complex Gaussian noise with zero mean 
and covariance matrix <j 2 Ik, where the scalar a 2 is the noise 
power. 

We collect N i.i.d observations from model (1) to a K x TV 
matrix X = [xi, . . . , Xjv], The problem of interest is to use 
the data matrix X to decide whether there are primary users. 
For ease of analysis we make the following assumptions 

1) The channel H is constant during sensing time. 

2) The primary user's signal follows an i.i.d zero mean 
Gaussian distribution and is uncorrected with the noise. 

In the absence of primary users, the sample covariance matrix 
R = XX^ follows an uncorrected (white) complex Wishart 
distribution Wk (N, X) with population covariance matrix 

S :=E[XX t ]/iV = cr 2 lK. (2) 

In the presence of primary users, by assumptions 1) and 
2), the sample covariance matrix R follows a correlated 
complex Wishart distribution. The correlation is induced by 
the presence of the signals, and it is characterized by the 
population covariance matrix 



-,K 



S = C7 2 I K +^ 7i h i ht, 



(3) 



where 7, := E[sjsJ] defines the transmission power of the i- 
th primary user. The received SNR of primary user % across 
the K sensors is SNR, := 7 i ||h i || 2 /<7 2 . Finally, we denote 
the ordered eigenvalues of the sample covariance matrix R by 

< \ K < ■ ■ ■ < Ai < 00. 

B. Test Statistics 

The differences between the population covariance ma- 
trices (2) and (3) can be explored to detect the primary 
users. This detection problem can be formulated as a binary 
hypothesis test, where hypothesis Ho denotes the absence of 
primary users and hypothesis Hi denotes the presence of 
primary users. Declaring wrongly Ho, or declaring correctly 
Hi, defines the false alarm probability Pf a , and the detection 
probability P^, respectively. Since the sample covariance ma- 
trix R is a Wishart matrix, it is a sufficient statistics for the 
population covariance matrix S [14]. This leads to various test 
statistics as functions of R with different assumptions on the 
number of primary users P, and the knowledge of the noise 
power a 2 . 

In the case of a single primary user, P = 1, the hypothesis 
test can be expressed as 



Ho 

Hi 



X = a 2 \ 
X = cr 2 I 



K 



K 



7ihih{. 



(4) 
(5) 



Assuming known noise power a 2 , the Largest Eigenvalue 
based (LE) detection XLe := Ai is shown to be optimal under 
the Generalized Likelihood Ratio Test (GLRT) criterion [2]. 
Here we concentrate on unknown noise power. Then the 
optimal detector in the GLRT sense is the Scaled Largest 



Eigenvalue based (SLE) detection T SLE :— Xi/^2 i=1 Aj. The 
SLE detector is first proposed in the context of spectrum 
sensing in [4] and analyzed in [5-8]. 

It can be verified that the matrix Ei=i 7jhjht in (3) is 
positive definite, i.e. Ei=i 7ihi n ! y 0. Considering this fact, 
in the presence of possibly multiple primary users, when P > 
1 but not known a priori, the hypothesis test is expressed as 



Ho 
Hi 



£ = o- 2 I 
X>~er 2 I 



(6) 

(7) 



where the noise power a 2 is assumed to be unknown. Essen- 
tially, we are testing a null hypothesis X = o 2 \k against all 
the other possible alternatives of X, i.e. the hypothesis test 
is blind to P. In the statistics literature, this hypothesis test 
is known as the sphericity test. The corresponding optimal 
detector under the GLRT criterion is the so-called Spherical 

Test based detection Tst := IIi=iAi/ (E»=iA»/-K') ■ I n 



the context of spectrum sensing, the ST detection is first 
proposed in [11] and analyzed in [12]. Although in general 
the ST detector achieves good performance, it is not the best 
one in the low SNR regime. A test statistics that is optimal 
in detecting small deviations from Ho is the so-called John's 
detection 



T,:= 



tr(R 2 



-,K 



X)»=i \ 



(*( R )) (eLa,)' 



(8) 



which is first considered by S. John [15]. A more rigorous 
derivation of the test statistics (8) can be found in [16], where 
the resulting test procedure is 



■Hi 

Tj ^ C, 

■Ho 



(9) 



C being a threshold. The optimality property of the Tj detector 
in detecting small derivations is known as the locally best 
invariant property. Mathematically, it means that for every a 2 
and for every other test T (say), there is a neighborhood of 
o 2 \k such that Tj achieves no worse performance than T 
does [16]. 

Besides the ST and John's detectors, another detector in the 
presence of multiple primary users is the Eigenvalue Ratio 
based (ER) detection T ER = Xi/X K [9, 10]. The ER detector 
is not constructed from any optimality considerations, thus its 
performance is substantially worse than those of the ST and 
John's detectors [6, 12]. Finally, we note that no eigenvalue 
decomposition is needed in implementing John's detector as 
opposed to most of other eigenvalue based detectors. 

III. Performance Analysis 

In this section we first derive an exact expression for 
the moment of Tj. Based on this result, we construct an 
approximation to the distribution of Tj, which leads to closed- 
form formulae for the false alarm probability and the decision 
threshold. 



A. Exact Moment Expression 

Under Ho, by following the similar argument for the real 
Wishart case [17] it can be easily verified that the random 

variable (Ei=i^») i s independent of the random variable 
of interest 



T, = 



(eL. 



e [l/tf, l] 



(10) 



K 



By this independence, the m-th moment of JZ i=1 A 2 equals 




E[Tj m ]E 



A' 



2m 



E^ 



A 



E^ 

vi=l 



(11) 



(12) 



^K 



The random variable 2tr(R) = 2^ i=1 Ai follows a Chi- 
square distribution with 2KN degrees of freedom. By using 
the moment expression for Chi-square distribution [18] (Eq. 
(2.35)), the 2m-th moment of 5Z i=1 A^ is obtained as 



E 



A 



E^ 



T(2m + KN) 
T{KN) 



(13) 



^K 



The next step is to calculate the moment of X^ili A 2 , which 
is given by the following result. 

Proposition 1. The m-th non-negative integer moment of the 
random variable ^2 i=1 A 2 equals 



E 



K 



E* 2 



E 



i\C 



oi!-- -a K \ 



K 

Yl {2 aj - 2a, h + j - i) J| r(2a 4 + N - K + i), (14) 

l<i<j<K j=l 

where the sum is over all the non-negative integer solu- 
tions of <X\ + • • • + a^ = rn and the constant C = 

The proof of Proposition 1 is in Appendix A. Inserting (14) 
and (13) into (12), the m-th moment of random variable Tj, 
denoted by M m , equals 

CT(KN) — m\ 



M m := 



E 



T(2m + KN) ^ ai\---a K \ 

K 



Yl {2a 3 - 2a t + j - i) J| r(2a 4 + N-K + i). (15) 

l<i<j<K j=l 

The sum over the partition a\ + ■ ■ ■ + sr- = 
to can be implemented by normal sums as 
Ql= o£a 2 =o •••Ea K - 1 =o wlth a K replaced 

by to — E;=i a i i n tne summand. Note that for real Wishart 
matrix, up to the second, fourth and sixth moment of Tj 



under Ho can t> e found in [17], [21] and [22], respectively. 
To the best of our knowledge, the derived moment expression 
of Tj (15) for complex Wishart matrix, which is valid for 
arbitrary non-negative moment, is new. 

B. Moment Based Approximation 

It is a standard technique in statistics to approximate 
some unknown distribution by a known one having the 
same support and moments. Motivated by the results for 
real Wishart case [22], in this work we choose a gener- 
alized Beta distribution with the same support as that of 
Tj to approximate the distribution of Tj. Specifically, the 
linear transform x = - — ~ ^ on a standard Beta density 2 
z a ~ x {\ — z)^~ 1 /B(a, (3), z £ [0,1] leads to a generalized 
Beta density 

Q-l 

('"[■'■- 4, \ (l-xf-\ (16) 



1 
K 



with the support x £ \}-/K, 1] an d the constant C" = 
K a+[)-i/ (B{a,0)(K - l)"^- 1 ). Since the m-th moment 
of a standard Beta random variable equals E[z m ] = 
{oc) m /(a + f3) m , where (a) m = T{pc + m)/T(a) defines 
the Pochhammer symbol, the m-th moment of the generalized 
Beta random variable is obtained by binomial expansion as 

'{K- l)z + V 



Eb" 



E 



— Y 

j=0 



i=0 



K 



(K - iyE[z* 



\i) (a + 0)i ' 



(17) 



(18) 



(19) 



where ( m ) denotes the binomial coefficient. In particular, the 
first two moments are 

(3) 2 + aK 2 



aK + f3 



{aK- 



P 



(20) 



[a + P)lC {a + /3){a + /3 + l)K 2 ' 

by matching them to the first two moments of Tj (15), the 
parameters a and j3 of the generalized Beta density (16) can 
be obtained 

(KM 1 -l)(KMi-KM 2 +M 1 -l) 



B 



(K-1)K(M 2 - 
{Mi - l)(KMi - KM* 



Mf) 



1) 



(21) 



(22) 



(K-1)(M 2 -M 2 ) 
As a result, the two-moment-based approximation to the CDF 
of Tj under Ho is 

Fj(y) « C f (x-iY {l-xf- 



'l/K 

B 



_1 dx (23) 



= 1 



K(l-y) . 



f3, a 



B(a,0) 

where y £ [1/K, oo) and B(x; a, b) 
is the lower incomplete Beta function. 



(24) 



/ X «°- 1 (1-«) 6_1 48 



2 B(a, 13) = r(a)F(/3)/r(a + 0) defines the Beta function. 



By (9), for a given threshold £ the two-moment-based 
approximation to the false alarm probability equals 



B 



PUC) = i - f,(0 



g(i-C) . 



/3,a 



(25) 



B(a,j9) 

To implement the proposed spectrum sensing algorithm, a 
decision threshold needs to be determined for a given detection 
requirement in a computationally affordable manner. Using 
the derived approximation to Pf a , an approximative decision 
threshold can be obtained by numerically inverting (25). Note 
that both a and f3 are elementary functions of the sensor size 
K and sample size N through (14). Moreover, if we further 
approximate a and /3 to their respective nearest integers, 
(25) reduces to a finite sum of polynomials in £, Thus the 
computational complexity of threshold calculation becomes 
quite affordable for on-line implementations. 

Here we note that the proposed two-moment-based Beta 
approximation corresponds to the simplest form of a general 
Jacobi polynomials based approximation. In the general frame- 
work, up to any n-th degree of Jacobi polynomials matching 
the corresponding n moments of Tj would be used. Since the 
random variable Tj is of a finite support, the Jacobi polyno- 
mials expansion for the distribution of Tj is exact according 
to Weierstrass approximation theorem [23]. Namely, when n 
goes to infinity the Jacobi polynomials based approximation 
represents the exact distribution of Tj. In practise, the choice of 
n reflects a trade-off between the approximation accuracy and 
the implementation complexity. In light of the good accuracy 
as shown in the next section, we consider n = 2 in this 
paper. The general rt-moment-based approximation, including 
the error analysis, can be easily obtained by following the 
procedures in [22, 24] . 

IV. Numerical Results 

In this section we first investigate the accuracy of the derived 
approximative false alarm probability by simulations. Then 
we compare the performance of John's detector with those of 
several detectors in realistic scenarios. The considered values 
of the parameters K and N in this section reflect practical 
spectrum sensing situations. The sample size N can be as 
large as a couple of hundred whereas the number of sensors 
K is at most eight due to physical constraints of the device 
size. 

A. Accuracy of the Approximative False Alarm Probability 

In Figure 1 we plot the approximative Pf a using (25) and 
the simulated Pf a as a function of the threshold. To quantita- 
tively show the approximation accuracy, we also calculate the 
average approximation error 3 of the proposed Pf a with respect 
to the exact Pf a as resulting from simulations. The results, 
summarized in the caption of Figure 1, show that the derived 
analytical Pf a matches the simulations well. 

3 Defined as (Y%=i \P&(Ci) " -JW(Ci)l) /n, where P fa denotes the ap- 
proximative false alarm probability, i.e. the RHS of (25), and n is the sampling 
size. In Figure 1 we assume uniform sampling in £ £ [0.125,0.3] with 
n = 10 7 . 
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Fig. 1. False alarm probability: analytical versus simulations. For K = 
8, N = 50, 100 and 200, the average approximation error on false alarm 
probability is respectively 5.94 X 10~ 8 , 5.03 X 10~ 8 and 5.07 X 10~ 8 . 
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Fig. 2. ROCs for P = 3. Subplot (a): SNRi = -6 dB, SNR 2 = -5 dB, 
SNR 3 = -4 dB with K = 4, N = 400; subplot (b): SNRi = 1 dB, 
SNR 2 = 2 dB, SNR 3 = 3 dB with K = 4, TV = 50. 



B. Detection Performance 

We compare the detection performance of John's detector 
with those of other known detectors by means of ROC curve. 
Since a ROC curve shows the achieved detection probability as 
a function of the false alarm probability, it reflects the overall 
detection performance for a given detector. Our focus here is 
detection in the presence of multiple primary users, thus we 
consider for comparison the ST detector. In addition the SLE 
detector, which is optimal for single primary user detection, 
is considered for comparison as well. Comparisons with the 
non-optimal ER detector and detectors that are sensitive to 
noise uncertainty [25], such as the LE detector and the energy 
detector [26], are excluded in this paper. For results in this 
direction, the readers are referred to [6,7, 12]. 

In Figure 2 we consider a scenario of three simultane- 



ously transmitting primary users with relatively low SNRs 
(SNRi = -6 dB, SNR 2 = -5 dB, SNR 3 = -4 dB using 
K = 4, N = 400) in subplot (a) and relatively high SNRs 
(SNRi = 1 dB, SNR 2 = 2 dB, SNR 3 = 3 dB using K = 4, 
N = 50) in subplot (b). Without loss of generality, we assume 
unit powers for the zero mean Gaussian signal and noise. 
The entries of the channel matrix H, which are fixed during 
sensing, are independently drawn from a standard complex 
Gaussian distribution. The same channel matrix is used in 
both subplots. The channel vector for each primary user is 
normalized as u^ = hj/||hj||. As a result, the population 
covariance matrix £ now equals £ — Ik + X)i=i SNRjUiU?. 
For the specific channel realizations considered in Figure 2, the 
eigenvalues 4 of £ in subplot (a) are [1.6225, 1.2217, 1.1213, 1] 
in subplot (b) are [4.0417, 2.2375, 1.56, 1], From Figure 2 (a) 
we observe that John's detector achieves the best detection 
performance in the low SNR case considered. However, when 
the SNRs increase we see from Figure 2 (b) that the ST 
detector performs better than John's detector. In both subplots, 
it is seen that the ST and John's detectors outperform the SLE 
detector. This is intuitively clear since the SLE detector is 
optimized for single primary user detection. Moreover, we see 
that as the SNRs increase the performance gap between the 
SLE detector and the multiple primary user detectors becomes 
larger, as expected. 

V. Conclusion 

In this paper, we investigated the sensing performance of 
John's detector, which a candidate detector in the presence of 
multiple primary users. John's detector is optimal in detecting 
small deviations of the covariance matrix from a matrix 
proportional to identity. Analytical formulae have been derived 
for the false alarm probability and decision threshold of 
John's detector. The derived results are simple to calculate and 
yield an almost exact fit to simulations. From the simulation 
setting considered, performance gain over several detection 
algorithms is observed in the low SNR regime. 

Characterizing 7} distribution under Hi, which leads to 
analytical results for detection probability, is the work in 
progress. 
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Appendix A 



K \2 



The non-negative integer moment of J2i=i A 

Under Ho, the joint density of the unordered eigenvalues 
Aj € [0, oo) for the sample covariance matrix R reads [19] 



C 



$(Ai,...,A„):=— |A(A)| 2 ;QA 



K 
2 TT \N-K -\ t 



(26) 



4 For the considered detectors, the test statistics depend on X only through 
the eigenvalues of S. 



where |A(A)| defines the determinant of Vandermonde matrix 
with i,j-th entry (A(A))j,j = \\~ , i,j = 1, . . . ,K and the 



tK 



constant C = [U4L1 r ( N ~ * + 1 ) r (- ft: - i + 1) 

Before we prove Proposition 1, we need the following two 
lemmas. 



K ,2o 



Lemma 1. The average value of the function Y\i = \\ at , 
where aiS are non-negative integers, equals 

I (frA 4 2a 'U(Ai,...,A K )dAi...dA K - 

J[o,°°) K Vti / 

J2 \Y{2a Vi +N-K + 1+J- 1)1^=1,...^ , (27) 



C_ 

Id. 



where v = v\ , . . . , vk defines a permutation of the integers 
1, . . . , K and the sum is over all the K\ permutations. 

Proof: Using the fact that a product of determinants 
equals the determinant of the matrix product and the fact that 
the determinant remains unchanged under transpose operation, 
we have 



|A(A)| 2 = |A(A)A'(A)| 



K 

E 
1=1 



a;+- 2 



(28) 



jj=i, 



.K 



By invoking the multi-linearity property of determinants, the 
Hankel determinant (28) above can be written as a sum of 
two determinants, where the first rows are [1, • • • , \ x ~ ] and 
[K — 1, • • • , X)i=2 \ _1 ] w i m th e res P ec t second to the last 
rows remain unchanged. By repeated use of the multi-linearity 
property, (28) can be written as sum of K K determinants, out 
of which K\ determinants with non-identical A index in rows 
give non-zero contribution, namely, 



K 

1=1 



E 



X i+j- 



,K 



(29) 



.K 



where v = v\,...,vk defines a permutation of the integers 
1,...,K and the sum is over all the K\ permutations. Inserting 
(29) into LHS of (27) and disregarding the constant C/K\ 
which will be re-installed, we have 



/ 

J[0,oo)' 



Z^ I "i \i,j=l,...,K 



K 



2a,i+N-K-\, 



dA,; 



E 






K 



[0, 



\ 2a Vi +N-K+i+j-2 _ Xv TT , x 

^Vi e * I Q-Ai 

\i,j=l,...,K J - J - 
i—1 



\T(2a^+N-K + i + j-l)\ 



*>j=ij 



,K ' 



where in the first equality we multipled each A 2ai+ c As 
with the row of |A' +:,_2 |. . , „ having the same A index 
and the second equality is achieved by first expanding the 
determinant using Leibniz formula, integrating term-wise and 
rewriting the integration results as a determinant. This com- 
pletes the proof. ■ 



Note that Lemma 1 can be considered as an extension to 
the Selberg type integral considered in [20] (Eq. (17.6.5) and 
Eq. (17.8.1)). 



Lemma 2. The following determinant can be simplified to 

K 

|r(&i + j-i)| iJ=Wf = " 



II (b,-h)l[r(b t ), (30) 

l<i<j<K j=l 



where hi is a positive integer. 

Proof: We first realize that from each row the term T(bi) 
can be factored out, namely, 

K 

\nb i +j-i)\ i j =1 ,...,K = \(pi) J -i\ iJ=1 ,... tK i[nb i ), 

(31) 
where (b)j = T[k=o(P+k). By extracting from the i-th column 
a suitable linear combination of previous i — 1 columns, the 
determinant 



1(6.) 



J- 1 \i,j=l,...,K 



br 1 



i,j=l,...,K 



(32) 



which is a Vandermonde determinant Y[i<i<j<K(bj~ &»)■ This 
completes the proof. ■ 

We are now in a position to prove Proposition 1. By using 
the multinomial expansion 



K 



E** a 



E 



ml 



K 



n*? 



and Lemma 1, we have 

( K \ " 



a\\ ■ ■ ■ «;<•! 

aiH \-a,K— Tn i—1 



(33) 



E 



E^ 2 



E 



i\ c 



ai _| [-aft:— m 



oi!---ajf!iin 



j; |r(2o Vi +JV-lif + * + j- l)| iji=1 «: ■ 



(34) 



For any given permutation u, it is observed that the sum over 
Oi + ■ ■ ■ + a-K = m is symmetric in the sense that one can 
arbitrarily permute the index of a without changing the value 
of this sum, in particular the following permutation of a holds, 



E 

oiH \-a,K- 


to! |r(2a„ s +N - K + i + j 


*-)\i,j=l,...,K 


ai! • • • a*-! 

-m 

to! |r(2a, + TV - K + i + j - 


*-)\i,j=l,...,K 



aiH ^aK— m 



di\ ■ ■ ■ Ojf! 



Since the number of possible permutations is K\, (34) now 
equals 



E 



m\C |r(2a, + N-K + i + j- 1)| . . x 



.K 



aiH \-aK— rn 



oi!-- 'Os-! 



(35) 

Using Lemma 2 with 6, = 2a^ + N — K + i completes the 
proof of Proposition 1 . 
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